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ABSTRACT 

Phenomenological relations exist between the peak luminosity and other observables 
of type la supernovae (SNe la), that allow one to standardize their peak luminosities. 
However, several issues are yet to be clarified: SNe la show color variations after the 
standardization. Also, individual SNe la can show residuals in their standardized peak 
absolute magnitude at the level of ~ 0.15 mag. In this paper, we explore how the color 
and luminosity residual are related to the wavelength shift of nebular emission lines 
observed at > 150 days after maximum light. A sample of 11 SNe la which likely suffer 
from little host extinction indicates a correlation (3.3(t) between the peak B — V color 
and the late-time emission-line shift. Furthermore, a nearly identical relation applies 
for a larger sample in which only three SNe with S — y > 0.2 mag are excluded. 
Following the interpretation that the late-time emission-line shift is a tracer of the 
viewing direction from which an off-centre explosion is observed, we suggest that 
the viewing direction is a dominant factor controlling the SN color and that a large 
part of the color variations is intrinsic, rather than due to the host extinction. We 
also investigate a relation between the peak luminosity residuals and the wavelength 
shift in nebular emission lines in a sample of 20 SNe. We thereby found a hint of a 
correlation (at ^ 1.6cr level). The confirmation of this will require a future sample 
of SNe with more accurate distance estimates. Radiation transfer simulations for a 
toy explosion model where different viewing angles cause the late-time emission-line 
shift are presented, predicting a strong correlation between the color and shift, and a 
weaker one for the luminosity residual. 

Key words: supernovae: general - cosmology: distance scale - cosmology: cosmolog- 
ical parameters 



* Partially based on observations obtained at the Gemini Obser- 
vatory, Cerro Pachon, Chile (Gemini Programs GS-2008B-Q-8, 
GS-2008B-Q-32, GS-2008B-Q-56, and GS-2009B-Q-40). 



1 INTRODUCTION 

Type la supernovae (SNe la) are used to measure cosmolog- 
ical parameters and study the nature of tlie dark energy (see 
Leibundgut 2008, and references therein). Thanks to the uni- 
formity of their peak luminosities, once a phenomenological 
relation between the light-curve shape and the peak lumi- 
nosity is applied (hereafter the light-curve correction, or the 
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Phillips relation), they can be accurately used as cosmo- 
logical standard candles (Phillips 1993; Hamuy et al. 1996; 
Phillips et al. 1999) . Their colors are also known to correlate 
with light-curve shape (Tripp 1998; Tripp & Branch 1999; 
Phillips et al. 1999). In addition to light-curve shape and 
color, several other observables, mostly related to spectral 
features, have been shown to correlate with the SN peak 
luminosity (e.g., Nugent et al. 1995; Mazzali et al. 1998; 
Benetti et al. 2005; Bongard et al. 2006; Hachinger et al. 
2006; Foley et al. 2008). 

Currently there are several issues yet to be clarified. One 
central issue is that the intrinsic color variations of SNe la 
have not yet been fully understood. After application of the 
relation between color and light curve shape, there remain 
variations in the color excess of SNe la. So far, it has been 
practically impossible to discriminate between the contribu- 
tions from a possible 'residual' intrinsic color, which does 
not correlate with the light curve shape, and that from the 
extinction within the host or the environment around the 
SN. This issue could be related to the fact that when the 
dispersion of the Hubble diagram is minimized with Rv be- 
ing treated as a free parameter, one obtains low values of Rv 
between 1 — 2. However, as shown by Folatelli et al. (2010; 
hereafter FIG), when one compares the colors or color ex- 
cesses of normal SNe la, a more typical Milky Way-like value 
of the reddening law is obtained, i.e. Rv ~ 3. FIO argued 
that this apparent discrepancy suggests that there is an in- 
trinsic color variation within SNe la that correlates with lu- 
minosity, but is independent of the light curve decline-rate 
Amis (B)I3 In the present study, we adopt an Rv = 1.72 as 
derived from FIO (i.e. Calibration 7 of Table 9). 

Moreover, after application of the light-curve correction, 
Hubble diagram residuals at the level of 0.15 mag exist 
for individual SNe la (e.g., Phillips et al. 1999; Prieto et al. 
2006; Jha et al. 2007; Hicken et al. 2009ab). This is one of the 
issues that presently limit the precision in using SNe la to 
constrain the value of the equation-of-state parameter of the 
dark energy (e.g., Hicken et al. 2009b; see also Wood-Vasey 
et al. 2007, Kessler et al. 2009 for the current status of the 
precision in SN la cosmology) . Several suggestions have been 
made for a secondary parameter that may provide a more 
accurate luminosity calibratioij3 (or on parameters already 
including the effect of the second parameter). Suggestions 
include (i) metallicity (Gallagher et al. 2005; Timmes et al. 
2003; Mazzah & Podsiadlowski 2006; Hoflich et al. 2010, 
but see also Howell et al. 2009; Neill et al. 2009; Yasuda & 
Fukugita 2010), (ii) high- velocity spectral features (Wang et 
al. 2009b), (iii) spectral flux ratios (Bailey et al. 2009; Yu et 
al. 2009), and (iv) the mass and/or the morphological type 
of the host galaxy (Kelly et al. 2010; Lampeitl et al. 2010; 
Sullivan et al. 2010). 

An interesting possibility for the origin of the diverse 
properties of SNe la was recently suggested by Kasen et al. 
(2009) theoretically and by Maeda et al. (2010ab, hereafter 
MlOa and MlOb) observationally, namely an asymmetry in 

^ Ami5(_B) is the magnitude difference in _B-band between max- 
imum brightness and 15 days later. 

^ Indeed, this is the 'third' parameter, since Hght curve fitting 
methods usually use two parameters, i.e., the light curve shape 
and color. In this paper, we simply call the additional parameter 
the 'second' parameter following the convention. 



the SN explosion combined with the observer viewing angle. 
In particular, MlOa identified potential signatures of asym- 
metry in a number of SNe la, based on the observed wave- 
length shift of late-time emission lines (see §2 for more de- 
tails). MlOab showed that the required configuration is qual- 
itatively consistent with the expectation from a deflagration- 
to-detonation transition scenaricjfl if the first thermonuclear 
sparks are ignited offset from the centre of the progenitor 
white dwarf. 

MlOb suggested that the viewing angle effect is a proba- 
ble origin of the spectral evolution diversity of SNe la. Differ- 
ent SNe la show different velocity gradients (wsi), defined as 
the speed of the decrease in the Si II absorption velocity after 
maximum brightness (Benetti et al. 2005; see also Branch et 
al. 1988) . SNe are divided into high- velocity-gradient (HVG) 
(■usi > 70 km s^^) and low- velocity-gradient (LVG) objects 
(■usi < 70 km s^^). MlOb argued that different velocity gra- 
dients are a consequence of different viewing directions from 
which the SN is observed. It has been indicated that LVG 
and HVG SNe may show different properties in their intrin- 
sic colors (e.g., Pignata et al. 2008) and that their luminosi- 
ties may have to be calibrated in a different manner (Wang 
et al. 2009b). Here we revisit this color issue in the context 
of our new interpretation of LVG and HVG SNe. 

In this paper, we explore whether the late-time emis- 
sion line shift, and thereby the observer viewing angle on an 
asymmetric explosion, is related to the intrinsic color and 
the luminosity residuals of SNe la after application of the 
Phillips relation. We find a correlation between the color at 
maximum brightness and the nebular emission line shift. We 
also investigate a possible relation between the luminosity 
residuals and the nebular line shifts, but since our sample is 
small the significance is not overwhelming. We then investi- 
gate the ramification of the viewing angle on the luminosity 
and color calibrations with the help of multi-dimensional 
radiation transfer calculations, and find that the predicted 
effect is qualitatively consistent with the trends seen in the 
data. 

The paper is organized as follows. In §2, we summarize 
the findings of MlOa regarding the asymmetry in SNe la, 
which are then used throughout the present paper. In §3, 
we present details of the sample of nearby SNe la consid- 
ered in this study. In §4, we discuss how the viewing angle 
is related to the intrinsic color of SNe la. In §5, we discuss 
the procedures to estimate the intrinsic absolute magnitude 
and subsequent residuals. In §6, we compare the late-time 
emission line shifts with the luminosity residuals. In §7, we 
investigate the effect of the viewing angle on the peak bright- 
ness and color by simulating light curves for kinematic off- 
centre toy models. In §8 the paper is closed with conclusions, 
discussion and future perspectives. 



2 ASYMMETRY IN TYPE lA SNE 

In this section, we summarize the findings presented in MlOa 
regarding asymmetry in SNe la. MlOa suggested that the 

In the deflagration-to-dctonation transition scenario, the ther- 
monuclear sparks first trigger the deflagration flames which travel 
subsonically, and then the flames subsequently turn into a super- 
sonic detonation flame (Khokhlov 1991). 
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Figure 1. [Fc II] A7155 and [Ni II] A7378 in late-time spectra of SNe la (see references in Table 1). The spectra have been redshift- 
corrected, and then converted to a velocity assuming 7378 A as the zero-velocity. The rest positions of [Fe II] and [Ni II] are denoted 
by dotted and dashed lines, respectively. The color coding indicates the early-phase 'velocity gradient' (see §1; Benetti et al. 2005): 
Low velocity gradient SNe (LVG SNe) in green, high velocity gradient SNe (HVG SNe) in magenta. SN 2007on is a fast declining SN 
1986G-like SN la (Morrell ct al. 2007) and is shown in black. SN 2006dd, without a measured velocity gradient, is shown in blue. SN 
2004dt, which was argued to be a peculiar outlier based on late-time spectra (MlOb), is shown in orange. 



innermost region of the ejecta, filled with stable Fe-peak 
elements, is generally offset, and that an observer's viewing 
angle can be traced by shifts in the central wavelengths of 
nebular emission lines of [Fe II] A7155 and [Ni II] A7378 at 
late phases. The argument is as follows: 

(i) Using a sample of late-time (i.e., at least 100 days after 
maximum brightness) SN la spectra, mostly drawn from the 
SUSPECT databaseQ MlOa found that in cases where the 
[Fe II] A7155 and [Ni II] A7378 lines were detected, they 
exhibited measurable shifts in their central wavelength with 
respect to the expected rest wavelength, as shown in Fig. 
1. There is no clear correlation between the phase of the 

^ The Online Supernova Spectrum Archive, SUSPECT, is found 
at |htt p : / /bruford . nhn . ou . edu/ ~ suspect / 1 . 



nebular spectrum and the shift of the [Fe II] A7155 and [Ni 
II] A7378 lines. In addition, no significant evolution is seen in 
the measured line shifts for individual SNe la at > 150 days. 
To illustrate these characteristics, we show the wavelength 
shifts (converted to a velocity) of these lines in Fig. 2. 

(ii) On the other hand, the strongest lines at late phases, 
i.e., the [Fe III] blend at ~ 4700 A and the [Fe II]/[Fe III] 
blend at ~ 5250 A, behave similarly for all SNe la. Since 
the [Fe III] blend at ~ 4700 A is stronger and its blended 
nature can be better handled (MlOa), we show its temporal 
evolution in Fig. 2. There is a clear correlation between the 
phase the nebular spectrum and the central wavelength of 
the [Fe III] blend at 4700 A, unlike for [Fe II] A7155 and [Ni 
II] A7378. The shifts of [Fe III] at 4700 A evolve from the 
blue to the rest wavelength for all SNe. 

(iii) The temporal behavior seen in Figure 2 suggests that 
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Figure 2. Observed wavelength shifts for (a) the [Fe II] A7155 and [Ni II] A7378 lines, and (b) the [Fe III] blend at 4700 A, plotted as 
a function of time since _B-band maximum. The wavelengths have been converted to velocities, assuming the expected rest wavelengths 
of the lines as the zero-velocity. The same objects are connected by lines. The data are from our present sample (see §3, Table 1 and 
references therein). 



the shifts of [Fe 11] A7155 and [Ni II] A7378 trace the line-of- 
sight velocity of the emitting material. As a consequence, the 
diversity in the wavelength shifts indicates that the distri- 
bution of the emitting material is offset with respect to the 
centre of the SN ejecta. Note that spherically symmetric ex- 
plosions with different expansion velocities should produce 
different line widths, but never a shift in the central wave- 
length. As a result of the offset, the observed diversity arises 
from different observer viewing angles. On the contrary, the 
observed features of the [Fe III] blend at 4700 A suggest that 
the [Fe III] blend traces the distribution of the emitting ma- 
terial only at sufficiently late epochs, and the small diversity 
indicates that the emitting material is distributed more or 
less spherically. 

MlOa argued that these characteristics can be natu- 
rally explained by a deflagration-to-detonation explosion 
scenario, if the distribution of the inner deflagration ash 
has an offset and the outer detonation ash is distributed 
in a roughly spherically symmetric way. They proposed that 
these lines are emitted from different regions - an outer, 
relatively low-density region dominated by radioactive ^®Ni 
(which decays into ^^Co and then into ^®Fe, powering the 
SN light curve), and an inner, relatively high-density region 
dominated by stable Fe-peak elements, i.e., ^*Ni, ^^Fe, and 
^''Fe (see also Mazzali et al. 2007). The former and latter cor- 
respond to the detonation and the deflagration ashes, respec- 
tively. For such a configuration, the outer region should be 
in a relatively high ionization state (i.e., doubly- ionized) and 
at high temperature (electron temperature Te > 10, 000 K). 
On the other hand, the inner region should be in a low 
ionization state (i.e., singly-ionized) and at low tempera- 
ture (Tc - 2000 - 7000 K). This stems from (a) the ion- 
ization balance and (b) the thermal balance. Namely, (a) 
ni+i/rii oc J^n~^, where rn is the density of the i-th ioniza- 
tion state, J-Y is the radioactive energy input from the decay 
chain ''"Ni ^'^Co ->■ ^^Fe, and ric is the electron density. 



(b) e""^™''"^" cx J~,n~^ng^ , where Tex is the excitation tem- 
perature of the line, and no is the population of the lower 
level. 

The [Fe III] blend at 4700 A is emitted by Fe++ with 
high Tox, and thus the outer region dominates the emission 
of this blend. This blend is thus attributed to emission from 
the detonation ash in the deflagration-to-detonation transi- 
tion scenario. The opposite is true for [Fe II] A7155 and [Ni 
II] A7378, and thus these lines are mostly emitted from the 
inner region, i.e., the deflagration ash. Combined with the 
phenomenologically derived distribution of material emit- 
ting these lines (see above) , this argues that the deflagration 
ash is located offset from the centre of the progenitor star, 
while the detonation ash is distributed spherically. 

Although these arguments are based on a few lines in 
the optical wavelength regime, this interpretation predicts 
that, at late phases, the lines from high ionization states 
and/or with high excitation temperature should show virtu- 
ally no shift, while those from low ionization states and/or 
with low excitation temperature show a wavelength shift de- 
pending on the viewing direction. MlOa showed that this is 
the case for SN 2003hv, which has a well observed spec- 
trum at late phases, covering spectral lines from the optical 
all the way to the mid-infrared (IR) (Motohara et al. 2006; 
Gerardy et al. 2007; Leloudas et al. 2009). The size of the 
offset derived for SN 2003hv was found to be ~ 3, 500 km s^^ 
(MlOa). Note that this behavior in various lines from opti- 
cal through mid-IR provides an additional argument against 
line blending and radiation transfer as a cause of the ob- 
served line shifts. 

This finding was interpreted as evidence that the initial 
deflagration flame proceeds in an asymmetric way, having 
an offset with respect to the explosion centre, in the context 
of a deflagration-to-detonation transition explosion scenario 
(Khokhlov 1991; Yamaoka et al. 1992; Woosley & Weaver 
1994; Iwamoto et al. 1999; Ropke & Niemeyer 2007; Kasen et 
al. 2009; Seitenzahl et al. 2010). Affhough the details of the 
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Table 1. SN la sample 



SN 


Host 


z 


(mag) 




(mag) 


(mag) 


B - V 
(mag) 


References^ 


"nob 
(km s^l) 


References^ 


1990N 


NGC 4639 


0.003395 


31.71 ± 0.15 


(KP) 


1.07 ± 0.05 


12.62 


0.037 


L98 


— 126 ± 600 


G96 


1994D 


NGC 4526 


0.001494 


30.98 ± 0.20 


(SBF) 


1.32 ± 0.05 


11.83 


-0.080 


R95, P96 


— 2220 ± 220 


G96 


1997bp 


NGC 4680 


0.008312 


32.68 ± 0.27 


(HF) 


0.97 ± 0.2 


13.78 


0.16 


R05 


2539 ± 410 


MlOb 


1998aq 


NGC 3982 


0.003699 


31.56 ± 0.08 


(KP) 


1.12 ± 0.05 


12.42 


-0.11 


R9g 


— 1106 ± 286 


BOS 


lyyobu 


INGC 33DO 


0.00z992 


OU. 11 ± U.zO 


(KP) 


1 .06 ± O.UO 


11.80 


0.34 


J 99, rlOO 


1 o An _1_ 1 *7 1 

— looy ±1/1 


COl 


2000CX 


NGC 524 


0.007935 


32.63 ± 0.27 


(HF) 


0.93 ± 0.04 


12.99 


0.10 


LOl 


-244 ± 600 


C03, S04 


2001cl 


NGC 1448 


0.003896 


31.23 ± 0.45 


(TF) 


1.13 ± 0.04 


12.69 


0.068 


K03 


993 ± 152 


M05 


2002bo 


NGC 3190 


0.004240 


31.70 ± 0.24 


(sbf£ 


1.16 ± 0.03 


13.50 


0.44 


B04 


2350 ± 100 


SOS 


2002dj 


NGC 5018 


0.009393 


32.82 ± 0.25 


(HF) 


1.08 ± 0.05 


13.85 


0.063 


P08 


2090 ± 423 


P08 


2002cr 


UGC 10743 


0.008569 


32.87 ± 0.25 


(HF) 


1.32 ± 0.03 


14.10 


0.16 


P04 


797 ± 600 


K05 


2003du 


UGC 9391 


0.006384 


32.42 ± 0.30 


(HF) 


1.02 ± 0.05 


13.54 


-0.089 


A05, S07 


-471 ± 265 


S07 


2003hv 


NGC 1201 


0.005624 


31.37 ± 0.30 


(SBF) 


1.61 ± 0.02 


12.49 


-0.035 


L09 


-2677 ± 457 


L09 


2004dt 


NGC 799 


0.019730 


34.55 ± 0.12 


(HF) 


1.21 ± 0.05 


15.25 


-0.053 


A07 


-1551 ± 600 


A07 


2004eo 


NGC 6928 


0.015701 


34.12 ± 0.14 


(HF) 


1.45 ± 0.04 


15.02 


0.064 


P07a 


-14 ± 600 


P07a 


2005cf 


MCG-1-39-3 


0.006461 


32.19 ± 0.33 


(HF) 


1.07 ± 0.03 


13.26 


-0.017 


P07b, WOO 


324 ± 600 


L07 


2006X 


NGC 4321 


0.005240 


30.91 ± 0.20 


(KP) 


1.31 ± 0.05 


13.96 


1.34 


W08, Y09 


1331 ± 164 


W08 


2006dd 


NGC 1316 


0.005871 


31.26 ± 0.10 


(PNLF) 


1.08 ± 0.03 


12.32 


-0.060 


SlOa 


-1569 ± 142 


SlOa 


2007OI1 


NGC 1404 


0.006494 


31.45 ± 0.19 


(SBF) 


1.55 ± 0.01 


12.98 


0.12 


SlOb 


272 ± 600 


SlOb 


2007sr 


NGC 4038 


0.005477 


31.51 ± 0.17 


(TRGB) 


1.12 ± 0.01 


12.54 


0.12 


SlOb 


1754 ± 198 


SlOb 


2009ab 


UGC 2998 


0.011171 


33.29 ± 0.20 


(HF) 


1.25 ± 0.02 


14.60 


0.058 


SlOb 


-634 ± 486 


SlOb 



Distance modulus from KP Cepheid measurements (KP), Surface Brightness Fluctuations (SBF), the ncar-IR TuUy-Fisher 
relation (TF), the planetary nebulae luminosity function (PNLF), the tip of the giant branch (TRGB), or the host-galaxy 
recession velocity corrected for the Virgo infall (HF). See §5 for details and references. 

^ Obtained from the literature, except for SNe 2007sr, 2007on, and 2009ab (Stritzingcr ct al. in prep.). For these 3 SNe, 
Ami5(B) estimates were obtained using the multi-color template light curve fitter SNooPy (Burns ct al. 2011). 
The values at maximum brightness. Corrected for Galactic extinction. 

References for the photometric properties. AOS. Anupama et al. (2005); A07. Altavilla ct al. (2007); BOS. Branch et al. 
(2003); B04. Benetti et al. (2004); COl. Cappellaro et al. (2001); C03. Candia et al. (2003); G96. Gomez et al. (1996); HOO. 
Hernandez et al. (2000); J99. Jha et al. (1999); K03. Krisciunas et al. (2003); K05. Kotak et al. (2005); L98. Lira et al. 
(1998); LOl. Li et al. (2001); L07. Leonard (2007); L09. Leloudas et al. (2009); M05. Mattila et al. (2005); MlOb. Maeda 
et al. (2010b); P96. Patat et al. (1996); P04. Pignata et al. (2004); P07a. Pastorello et al. (2007a); P07b. Pastorello et 
al. (2007b); P08. Pignata et al. (2008); R95. Richmond et al. (1995); R99. Riess et al. (1999); R05. Reindl et al. (2005); 
S04. SoUerman ct al. (2004); SOS. Stchle et al. (2005); S07. Stanishev et al. (2007); SlOa. Stritzingcr ct al. (2010); SlOb. 
Stritzinger et al. (in prep.); W08. Wang ct al. (2008); W09. Wang et al. (2009a); Y09. Yamanaka ct al. (2009) 
^ References for the late-time spectra. 

/ An SBF distance to NGC 3226, a member of the same group. 



ignition process are not yet fully understood, theoretically 
an ofF-centre ignition may be a natural consequence of con- 
vection within a progenitor white dwarf (Kuhlen et al. 2006). 
Maeda et al. (2010c) (hereafter MlOc) argued, based on 
their hydrodynamic and nucleosynthesis simulations, that 
a deflagration-to-detonation transition model with the ini- 
tial sparks ignited at an offset from centre will result in a 
configuration qualitatively consistent with the above find- 
ings. 



According to this interpretation and an extensive search 
of emission lines from the optical to mid-IR wavelengths, 
MlOa suggested that the best lines to use in the optical for 
probing the asymmetry in the innermost SN ejecta are [Fe 
If] A7155 and [Ni II] A7378 (as well as [Fe II] A8621; Leloudas 
et al. 2009). These are the only ones which satisfy the criteria 
that (a) they reflect the low ionization and low temperature, 
(b) they suffer from little blending, and (c) they are covered 
by most optical spectra. We therefore adopt these lines as 
diagnostics of the ejecta asymmetry and the viewing angle 
in the present study. 



3 SUPERNOVA SAMPLE 

In this paper, we investigate how an asymmetric explosion 
and different viewing angles influence the color and lumi- 
nosity of SNe la, by examining correlations between these 
quantities and the emission line shift in the late-time spec- 
tra. The sample of SNe la is thus limited by the requirement 
that late-time nebular spectra are available. Our initial late- 
time spectral data set comprised the same 20 SNe la used by 
MlOa, mostly drawn from the SUSPECT database, supple- 
mented by other late-time nebular spectra accessible to the 
authors. To include these spectra in the present analysis, we 
imposed some additional criteria: (i) the spectrum was taken 
at > 150 days after maximum brightness (see Fig. 2), (ii) 
eithM [Fe II] A7155 or [Ni II] A7378 could be identified, (iii) 
the spectra had to be of high enough quality that the central 
wavelength of either of these lines could be measured, (iv) 
the SNe were required to have good early phase B- and V- 
band photometry so that their light-curve parameters and 
color at maximum could be estimated, and (v) Ami5(B) 
was in the range between 0.7 to 1.7 mag in order to apply 
the well-developed relations between Ami5(_B) and the peak 
absolute magnitude/color (e.g., Phillips et al. 1999; FIO). 

From the initial data set, 13 objects fulfilled these crite- 
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ria. In addition, we included 7 SNe which satisfied our crite- 
ria, i.e. SNe 1997bp, 2002bo, 2006X, 2006dd, 2007on, 2007sr, 
and 2009ab. In 13 out of these 20 objects, both [Ni II] A7378 
and [Fe II] A7155 were securely identified (see below). In 
the other 7 events, only one of these lines was discernable. 
These 7 SNe la were included in our analysis because a mea- 
surement of their nebular line velocity shift was possible, 
although less secure than when both lines were identified. 
The sample of 20 SNe la used in the present study is listed 
in Table 1. Their late-time spectra in the wavelength range 
covering the [Fe II] and [Ni II] are shown in Figure 1. Details 
on SN 2006dd are presented in Stritzinger et al. (2010), while 
SNe 2007on, 2007sr, and 2009ab will be presented elsewhere. 

MlOa measured the wavelength shifts in [Ni II] A7378 at 
late phases for the majority of the SNe la listed in Table 1. 
In MlOb, they also measured the shift in [Fe II] A7155. In 
this study, we remeasured these line shifts and updated the 
velocity shift (uneb), as follows: when both [Ni II] A7378 and 
[Fe II] A7155 were identified, we fitted their central wave- 
lengths with Gaussian profiles simultaneously using the IRAF 
dehlend commandlf] The v^^h is taken to be the mean value 
of them (see also M10b)l3 As a default estimate, the con- 
tinuum is taken to be constant across the wavelength range. 
We estimated the error in the fit by varying the continuum. 
We further considered the errors arising from the difference 
in the velocity shifts of the two lines. The regions emitting 
these two lines are not exactly the same since the region 
emitting [Ni II] A7378 is attributed to be a product of the 
deflagration at the very beginning, while that emitting [Fe 
II] A7155 is a product from a later phase, but still before the 
ignition of the detonation (MlOa; see also MlOc for an ex- 
ample of a detailed nucleosynthesis study). Therefore, there 
could be intrinsic differences in the velocity shifts of these 
two lines. The final error bars are taken to be the larger one 
between the errors associated with the Gaussian fit and the 
difference in the two lines. 

In the majority of SNe, The shifts seen in [Ni II] A7378 
agree with those measured in [Fe II] A7155 to within an er- 
ror of at most ~ 600 km s~^. If only one of these lines could 
be identified, we estimated the velocity shift only from the 
single line using the (single-line) Gaussian fitting command 
within IRAF, and assumed that the error is 600 km s"\ In 
the case of the fast-declining 1986G-like SN 2007on (Morrell 
et al. 2007) and the peculiar (in terms of the late-time spec- 
tra) SN 2004dt (MlOb), we did not use the [Ni II] feature as 
it might be contaminated or originate from [Ca II] AA7291, 
7324 (Filippenko et al. 1992; Turatto et al. 1996). Indeed, we 
suspected that [Ca II] AA7291, 7324 might contribute much 
to the "[Ni II] feature" for SNe 1994D, 2003hv, 2004dt, and 
2007on, since the measured position of the [Ni II] feature 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

^ The two lines were not simultaneously fitted in the previous 
measurements, and this changed the exact values of J^nob- How- 
ever, the old and new measurements are mostly consistent within 
the associated errors. The only exception is SN 1990N, for which 
the host galaxy redshift was incorrectly treated in the previous 
measurement. We confirmed that the conclusions of MlOab are 
not affected by the difference in the measurement in iinob- 



0.2 



E 



-0.2 



- 


■ 1 ■ 1 ■ 






— • 1 


O 

F10 






■ 









0.8 1.0 1.2 1.4 1.6 1.8 
Am,j(B) [mag] 



Figure 3. The -Bmax — Vmax color of 11 "low-extinction" SNe as 
a function of Ami^{B) . The color is corrected for Galactic extinc- 
tion only. SNe are shown by different symbols defined as follows: 
Low velocity gradient SNe (LVG SNe) by green fiUod-squares, 
high velocity gradient SNe (HVG SNe) by magenta filled-circles, 
the fast-declining SN 2007on by a black open-circle. SN 2006dd, 
for which no information on its velocity gradient is available, is 
shown by a blue open-triangle. The relation between color and 
Ar?ii5(_B) as derived by FIO is shown by a dashed line. 

is close to the rest wavelength of the [Ca II]. Among these 
SNe, we regard that the [Ni II] identification is the case for 
SNe 1994D and 2003hv, since its velocity shift agrees with 
that measured from the [Fe II]. On the other hand, these ve- 
locities are inconsistent for SNe 2004dt and 2007on, thus we 
believe that [Ca II] is a more likely identification for these 
two SNe. This indicates that the strong [Ca II] in the late 
phase may be a property shared by a part of fast declining 
SNe and related objects, but not by Branch-normal SNe. 
When spectra at several epochs were available, we adopted 
v-nch estimated from the latest epoch for which the quality 
of the spectrum was sufficiently high, since a more secure 
estimate of the intrinsic velocity shift for the nebular emis- 
sion lines is possible using later phase data (§2). As shown 
in Fig. 2, v^eh measured for the same SN at different epochs 
agrees within the errors, at least for > 150 days. 



4 INTRINSIC COLOR VARIATIONS 

We first examine whether the intrinsic colors of SNe la are 
related to v-nch and therefore to the viewing direction of an 
observer in the asymmetric explosion scenario. For this pur- 
pose, we define the pseudo color (iJmax — V^ax) as the dif- 
ference in peak B and V magnitudes obtained at the time 
of maximum light in each bandpass. 

We first identify a "low-extinction" sample of SNe for 
which the extinction within the host galaxy is likely insignif- 
icant. We regard SNe as low-extinguished if (i) the host is 
an elliptical or SO galaxy or (ii) the SN was located in the 
outskirts, further from the centre than half of the apparent 
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Figure 4. The Bmax — Vmax color, after correction of the Ami5 (B) - color relation of FIO (c(Ar?ii5 {B)): see the main text), vs. Dncb- The 
dashed line is the best linear fit to the 11 low-extinction SNe. (a) shows the low-extinction sample, and (b) the entire sample (zoomed 
in the same region as in the panel (a)). The inset in (b) shows the entire sample, including also the 3 highly-reddened SNe la. SNe are 
shown by different symbols defined as follows: Low velocity gradient SNe (LVG SNe) by green fiUed-squares, high velocity gradient SNe 
(HVG SNc) by magenta fiUed-circles, and the fast-declining SN 2007on by a black open-circle. SN 2004dt, which was suggested to be a 
peculiar outlier based on its late-time spectrum, is shown by an orange open-square. SN 2006dd, for which no information on its velocity 
gradient is available, is shown by a blue open-triangle. 



radius of the host (defined as the length of the projected ma- 
jor axis at the isophotal level 25 mag arcsec"^; Paturel et 
al. 1991) as obtained from the Lyon-Meudon Extragalactic 
Database (Paturel et al. 2003). We also remove SN 2000cx, 
since this was reported to show pecuUar colors (Li et al. 2001; 
Candia et al. 2003). Our adopted "low-extinction" sample 
then comprises 11 objects. 

It has been shown that the colors of SNe la at maxi- 
mum brightness are related to Amis {B) such that SNe with 
a larger decline rate are intrinsically redder (e.g., Phillips et 
al. 1999). Figure 3 shows a comparison between Ami5(_B) 
and the pseudo color (Bmax — Kiiax) for our low-reddening 
sample, corrected for Galactic extinction adopting the val- 
ues derived from the dust maps of Schlegel et al. (1998) and 
Rv = 3.1. Following the suggestion by MlOb that the view- 
ing direction is related to the velocity gradient, we hereafter 
plot the SNe with different symbols depending on whether 
they belong to the LVG or the HVG groups. Although we do 
not see a strong correlation, our data are overall consistent 
with the Ami£,{B) - color relation of the "low-reddened" 
sample of FIO defined as c{Amis{B)) = (Bmax — Kiax)o = 
-0.016 + 0.12[Ami5(B) - 1.1] mag. 

Figure 4 shows a comparison between Dnob and pseudo 
color corrected for Galactic extinction and the color vs. 
Ami5(_B) relation of FIO. In other words, we are investi- 
gating a relation between Vneb and the 'color residual' af- 
ter correcting for the Amis{B) term. Figure 4 reveals that 
there is a clear correlation between v^cb and this color for 
the low-extinction sample. A linear fit to the 11 SNe re- 
sults in a chance probability of only 4.4 x 10"^ (i.e., 3.3cr 
significance). The standard deviation is 0.05 mag. Figure 4 
further shows that almost the same relation applies even 
when not restricted to the low-extinction SNe (but omitting 



3 evidently highly-reddened SNe). This suggests that this 
'^ncb - color relation is at least as important as the Ami5(_B) 
correction in the color. It further indicates that a large part 
of the color variations of SNe la can be attributed to in- 
trinsic variations due to different viewing directions, rather 
than the host extinction. The intrinsic color we derive can 
be expressed by the following equation: 

max — 'Knax)() — 0.016 -f0.12[Ami5(B) - 1.1] 

4- 0.047(Dncb/1000 km s'^) mag. (1) 

It has been suggested that some HVG SNe may be in- 
trinsically redder than LVG SNe (e.g., Pignata et al. 2008). 
Wang et al. (2009b) showed that SNe with higher Si II 
absorption velocity (which is a typical signature of HVG 
SNe) are intrinsically redder than those with lower velocity 
(roughly corresponding to LVG SNe) by ~ 0.1 mag. There- 
fore, we qualitatively expect that SNe with larger (positive) 
Vncb are intrinsically redder according to MlOb, which is 
consistent with the data shown here. This also explains why 
Wang et al. (2009b) obtained a smaller dispersion in the cal- 
ibrated luminosities by assuming a smaller Rv for SNe with 
higher velocity in the Si II absorption (i.e., roughly equiva- 
lent to HVGs). If these SNe are intrinsically redder than the 
others, then a smaller Rv mimics the effect of the viewing 
direction. 

Figure 5 shows a comparison between the velocity gra- 
dient and -Bmax — Vmax (corrected for Galactic extinction and 
the Ami5(_B) term). As suggested previously, HVG SNe are 
redder than LVG SNe. The Unear fit results in P = 0.0073 
(2.4(7) for the low-reddened SNe. We note that the distribu- 
tion can actually be better expressed by a bimodal distri- 
bution (HVG SNe and LVG SNe) rather than a continuous 
distribution. Velocity gradients are similar for all LVG SNe, 
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Figure 5. The Bmax — Vmax color, after the correction of the Ami5(i?) - color relation of FIO (c(Ami5(B)); see the main text), as a 
function of the velocity gradient i)gi. The dashed line is the best linear fit to the 10 low-extinction SNe (excluding SN 2006dd for which 
the velocity gradient is not available). The symbols and color coding indicate the velocity gradient near maximum brightness (sec Fig. 4 
caption) . 



whereas we can see a color variation as a function of iincb also 
within this group. The correlation with the velocity gradient 
is indeed weaker than the relation using Uncb- This suggests 
that Uncb is a better indicator of the intrinsic color than the 
velocity gradient. This mostly stems from the saturation of 
the velocity gradients for LVG SNe, which makes it difficult 
to derive the detailed viewing direction directly from the 
velocity gradient (MlOb). 

Our finding that color is strongly correlated with Uncb 
suggests that the viewing direction is indeed an important 
property which controls the color as theoretically investi- 
gated by Kasen et al. (2009) (see also Foley & Kasen 2010). 
MlOb attributed the variation in velocity gradients in dif- 
ferent radial density distributions at different directions; 
they suggested that the photospheric velocity is smaller and 
evolves more slowly for a direction closer to the offset igni- 
tion direction. A smaller photospheric velocity (uph) leads to 
a larger photospheric temperature (Tes) since Teff oc v^^ 
The larger Teff results in a bluer color (c.f., figs. 1 and 3 
of Nugent et al. 1995). Also, the difference in Tcfi has been 
suggested to be a main origin of different spectral features 
seen in HVG and LVG SNe (Hachinger et al. 2008; Tanaka et 
al. 2008). This qualitatively explains the tendency we have 
identified in the data. 

This result opens up a possibility to estimate the host 
extinction by including the effect of viewing angle on differ- 
ent SNe: we compute the color excess due to host extinction 
as E{B - V) = (Bmax - Kiax) - (Smax - Vmax)o, where the 
SN intrinsic color (Bmax — Vinax)o is expressed by Equation 
1. Note that some SNe then show small negative color ex- 
cesses, since Equation 1 expresses the mean behavior of the 
low-extinction SNe with the standard deviation of ^ 0.05 
mag. We have used this method to estimate the extinction 
within the host or the environment around the SN and have 
tabulated our results in Table 2. 

We then compare the derived color excesses to those 



previously estimated in the literature. We focus on a com- 
parison with the E{B — V) computed by Wang et al. (2009b), 
since these were derived in a homogeneous and systematic 
way (and there is a sufficient overlap with our samples), but 
in Appendix A we also compare our values with other es- 
timates of E{B — V) found in the literature (reaching the 
same qualitative conclusion). Figure 6 shows a comparison 
between E{B — V) as derived in this study and the values 
from Wang et al. (2009b), for 16 SNe which are common in 
the two studies. Figure 6 shows that the color excess we de- 
rived tends to be smaller, since a large part of the color vari- 
ations are interpreted to reflect the intrinsic color variations 
rather than the host extinction. Indeed, the slope between 
the two estimates is consistent with ^ 3/1.85 for objects 
with E{B — V) < 0.2, where 1.85 is the value of Ry obtained 
by Wang et al. (2009b) for their whole sample. This means 
that even if we assume Ry ~ 3 to convert our E{B — V) 
estimate to A\r, the extinctions we derive are mostly consis- 
tent with those derived by the other method. Exceptions are 
the heavily reddened SNe, especially SN 2006X. This means 
that at least part of the reason why a small Rv is preferred 
for the luminosity standardization might be attributed to an 
overestimate of the host extinction (FIO), although at least 
SN 2006X clearly requires small Rv- 



5 COMPUTING INTRINSIC LUMINOSITIES 
AND RESIDUALS OF TYPE lA SNE 

To investigate any residuals in the intrinsic luminosities of 
SNe la, it is crucial that the distance to each event is esti- 
mated carefully. The distance should be provided indepen- 
dently of the SN properties. Most of the SNe la in our sample 
are nearby (that is why late-time spectral observations were 
possible for these objects), and therefore the most accurate 
distance estimates to their host galaxies come from direct 
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Figure 6. Comparison of the extinction. The horizontal axis de- 
notes the values of E(B — V) derived in this study corrected for 
Vjieb (Table 2), while the vertical axis provides values from Wang 
et al. (2009b). The comparison is shown for 16 SNe which are 
common between the two studies. The symbols and color coding 
indicate the velocity gradient near maximum brightness (see Fig. 
4 caption). 



distance measurements such as Cepheid variables or the Sur- 
face Brightness Fluctuation (SBF) method. Although the 
Cepheid measurements, available for 4 SNe in our list, are re- 
garded as the most reliable estimates, there is disagreement 
in the different analyses of the Cepheid data. The Cepheid 
distance measurements have mainly been contributed by 
the Saha-Tammann-Sandage SN la Hubble Space Telescope 
Calibration Program (hereafter STS: Saha et al. 2006) and 
by the HST key project (hereafter KP: Freedman et al. 2001; 
Stetson & Gibson 2001). The KP distances are typically 
shorter by 0.2 — 0.3 mag than those of STS (see Riess et 
al. 2005; Wang et al. 2006, for a discussion). Unfortunately 
such a systematic difference is important for our study, since 
it is exactly the order of magnitude of the effect we hope to 
investigate. We have chosen to adopt the KP distances since 
they have been suggested to give smaller dispersion in the 
residuals (e.g., Wang et al. 2006), and since this is the dis- 
tance scale the SBF distances have been calibrated to (see 
below). In Appendix B, we discuss how our results would be 
modified if we instead had adopted the STS distances. 

In the case of the SBF measurements, we used the values 
given by Tonry et al. (2001), and included the 0.16 mag 
correction suggested by Jensen et al. (2003) to fit to the 
KP Cepheid zero-point. For the SN host galaxies not listed 
in Tonry et al. (2001), we used values given by Ajhar et 
al. (2001) in which the zero-point is calibrated in a similar 
manner. 

There are several objects which require additional con- 
siderations about the distance. 

• For SN 2000CX in NGC 524, an SBF distance of /x = 
31.74 ± 0.20 mag is available. However, this value is signif- 
icantly smaller than the Hubble fiow distance (see below) 



of 32.63 ± 0.27 mag. We suspect that the SBF distance to 
NGC 524 is not correct, and have decided to adopt the Hub- 
ble fiow distance for this SN (see also Appendix B). 

• Neither a Cepheid nor an SBF distance is available for 
SN 2001el. We therefore used a distance based on the near- 
IR TuUy-Fisher relation (Wilhck et al. 1997). 

• For SN 2002bo, direct distance measurements are not 
available. However, there are SBF distances measured for 
two possible members of the group (see Krisciunas et al. 
2004): NGC 3193 {^l = 32.50 ± 0.18 mag) and NGC 3226 
(31.70 ± 0.24 mag), which however, do not agree with each 
other. The latter value is consistent with the Hubble flow 
distance (31.44 ± 0.48 mag) and we have therefore decided 
to adopt this latter SBF value. 

• For SN 2006dd, Feldmeier et al. (2007) reported a dis- 
tance to NGC 1316 based on the planetary nebula lumi- 
nosity function (PNLF) as = 31.26 ± 0.10 mag, which 
Stritzinger et al. (2010) argued is more reliable than the 
SBF distance (31.50 ± 0.17 mag). Therefore, we adopt the 
PNLF distance in this paper. 

• In the case of SN 2007sr, Schweizer et al. (2008) rean- 
alyzed Hubble Space Telescope data used by Saviane et al. 
(2008) and obtained a TRGB (tip of the red giant branch) 
distance of /x = 31.51 ± 0.17 mag, which is in agreement 
with the Hubble fiow distance. We adopt this value as the 
distance to SN 2007sr. 

Finally, there are a number of events for which we had to 
rely on the Hubble fiow (HF) distance scale, since estimates 
based on the aforementioned methods are not available. 
Recession velocities were corrected for Virgo infall as ob- 
tained from the Lyon-Meudon Extragalactic Database (Pa- 
turel et al. 2003), and converted to a distance modulus using 
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Ho = 72 km s"^ Mpc"^ (Freedman et al. 2001). The correc- 
tion with respect to the heliocentric radial velocity due to 
the Virgo infall is typically less than 100 km s~^. In these 
calculations we assumed an uncertainty of 300 km s~^ to 
account for peculiar motions. 

Reddening due to dust is the other main problem to es- 
timate the SN luminosity. Although we suggest that a large 
part of the color variations is intrinsic at least if we omit 
highly-reddened SNe (§4), the corrected color excesses are 
not always negligible (Tab. 2). For Galactic extinction we 
adopt the dust maps of Schlegel et al. (1998) with Rv = 3.1. 
We have estimated the host extinction using the method de- 
veloped and described in §4 (see Equation 1). The associated 
E{B — V) color excess (Table 2) was then converted to Ay 
using the relation Av ~ RvE{B — V). The associated error 
is: 

a^,, = RWe(b-v) + [E{B - V)fol^ . (2) 

Estimating the uncertainty in Ry is not trivial. We have 
assumed the typical value derived by FIO, i.e. aa^ = 0.3. 
Following the procedure of the previous section, we adopted 
an error of 0.05 mag for E{B — V). We assumed Rv = 1.72 
since to derive the residuals we used the light curve correc- 
tion from FIO, who obtained this value as the best fitting 
parameter. In Appendix B, the effect of the uncertainty in 
Rv upon our results is further investigated. 

We note that there are two approaches in treating the 
color and extinction for SNe la (see, e.g., FIO for further 
discussion). In one approach, the color and the light curve 
shape are simultaneously varied to obtain the standardized 
luminosity (e.g., Guy et al. 2007). In the other approach, 
first the intrinsic color is associated with other observed pa- 
rameter(s) (e.g., Ami5(_B)) and any excess beyond this color 
is assumed to be caused by the host extinction (e.g., Phillips 
et al. 1999). We follow the latter approach in this paper. We 
however emphasize that our treatment takes into account 
the intrinsic color variations independent from Ami5(_B), as 
is done in the former, simultaneous-fit approach. 

The observed l/-band peak magnitude was converted to 
absolute peak magnitude [My) using the distance modulus 
and the extinction as described above. We have not included 
K-corrections, since the correction is at most ~ 0.01 mag at 
maximum brightness for SNe at the low redshifts consid- 
ered here (Hamuy et al. 1993; Nugent et al. 2002). For the 
residual, the difference between the absolute peak V mag- 
nitude {Mv) and the standardized magnitude predicted by 
the Phillips relation [Mv(Ami5(_B))] was computed. In ob- 
taining the "standardized" peak magnitude, we applied an 
updated relation presented by FIO (see their Table 9.). Here- 
after, we denote the residual by dm = Mv —Mv{/S.mis{B)). 
Table 3 summarizes the Mv values that we computed, along 
with the standardized peak magnitudes. Also shown are the 
estimated errors in the calibration. 

In the resulting distribution of My, the dispersion is 
at the level of ~ 0.28 mag. This value is larger than that 
typically derived from samples of SNe la (~ 0.1 — 0.2 mag; 
e.g., FIO). We note, however, that such studies make use of 
uniform data sets and minimize the dispersion by parame- 
ter fitting the entire data set, including simultaneous fit to 
light-curve shape (e.g., Ami5{B)) and color (i.e., Rv and 
intrinsic SN color). Given that we have compiled data from 
various sources and adopted an external relation (FIO) de- 



Table 3. Peak Luminosity Calibration 



SN Afv{Ami5(B)) My dm <TAy2. '^My^ 

(mag) (mag) (mag) (mag) (mag) 



1990N 

1994D 

1997bp 

1998aq 

1998bu 

2000CX 

2OOI0I 

2002bo 

2002dj 

2002er 

2003du 

2003hv 

2004dt 

2004eo 

2005cf 

2006X 

2006dd 

2007on 

2007sr 

2009ab 



Error in the extinction estimate. 

Total error including the extinction and distance uncertainties. 




-3000 -2000 -1000 1000 2000 3000 



Figure 7. My as a function of iincb- The symbols and color 
coding indicate the velocity gradient near maximum brightness 
(see Fig. 4 caption). 

rived from an independent data set, a larger dispersion is 
expected. More importantly, most of the SNe la in our anal- 
ysis are nearby and are not in the smooth Hubble flow, so 
that the uncertainty in the distance estimates is larger (due 
to peculiar velocities; see for example, fig. 19 of FIO). 



6 RELATION BETWEEN THE LATE-TIME 
EMISSION LINE SHIFT AND THE 
LUMINOSITY RESIDUALS 

In this section we explore whether there are any relations be- 
tween Vnah and the luminosity of SNe la at maximum bright- 
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Figure 8. The residuals dm as a function of ■Uneb (see Figure 4 
caption). Tlie black dashed line shows the best linear fit to the 
entire data set (Tab. 5). 



Figure 9. Fractional number distribution of the 20 SNe la as 
a function of dm. SNe la are divided into two categories: those 
showing a blueshift in the late-time emission lines (gray) and 
those showing a redshift (white). 



ness. We have first investigated whether Vnch correlates with 
the absolute magnitude Mv- In Fig. 7 it is shown that there 
is no apparent correlation between these two parameters. 
The observed distribution of Mv as a function of Vneb can 
emerge from a non-correlation at the 14% level (P = 0.14), 
as derived by our Monte-Carlo simulations (see Appendix 
C). There is a clear outlier in the plot, that is SN 2000cx 
(see also Li et al. 2001). Even if we omit this single SN from 
the fit, the fitting result is not improved (P = 0.19). 

In the next step, we investigated if v^eb correlates with 
the luminosity residual dm (Fig. 8). According to the Monte- 
Carlo simulations (see Appendix C), these two quantities are 
correlated at 1.6a level. The best linear fit to the 20 data 
points, expressed by dm = a (uncb/1000 km s~^) + /?, is 
given in Table 4. We note that this investigation is limited 
by the relatively large errors in dm associated with the dis- 
tance uncertainties to the SNe (§5 and Appendix B). If we 
arbitrarily choose to ignore these errors (that were computed 
with a range of different methods) and consider only the er- 
rors associated with the extinction (that were derived in a 
systematic way), we observe that the significance of the pro- 
posed correlation increases (2.4(7; P — 0.0092). Although the 
significance for the observed correlation is low, it deserves 
further study when a larger sample is available. 

There may be a concern that the apparent relation in 
the Vneb — dm could be dominated by single points, rather 
than expressing a general trend. To examine this, we ar- 
tificially removed single points from the data set, and de- 
termined how much the quality of the fit was affected. We 
thereby obtained P = [0.038 - 0.086] (1.4cr - 1.8cr), depend- 
ing on which of the 20 points was removed. Also, the result- 
ing fits were all consistent within the la error. 

We also split the SN sample into halves according to 
the value of Ami5(B), and repeated the fitting for both 
subsamples (Table 4). We find that the 'faint' group (with 



Ami5(B) ^ 1.13 ma^ shows a correlation with a signifi- 
cance similar to that using the entire sample. On the other 
hand, the correlation is not evident for 'bright' SNe with 
Ami5(B) ^ 1.12. We note that the best fit slope for the 
faint sample is steeper than the one for the entire sample, 
as expected from the behavior of the bright sample. A pos- 
sible interpretation is presented in §7. 

We have also investigated a tendency by dividing the 
sample depending on Wnob- Figure 9 shows the number frac- 
tion of SNe la as a function of dm. In this figure, the sub- 
sample with Vnch < (blueshift) and that with Vneb > 
(redshift) are marked differently. Fig. 9 shows that there is 
a tendency that SNe la having blueshifts/redshifts in the 
late-time emission lines are brighter/fainter than expected 
from the Phillips relation. 

The difference in the average values for dm in the 
'blueshift' group (A'^ — 11) and 'redshift' group (A'^ — 9) 
is ~ 0.25 mag. To examine the statistical significance, we 
performed Monte-Carlo simulations as follows. First, we ran- 
domly selected 11 out of the 20 SNe ("A group"), and then 
the remaining 9 SNe were labeled as "B group". We then 
calculated the difference in the average values of dm in the 
two groups. With 10^ Monte-Carlo realizations, we find that 
the probability, P, that dm (A) — dm (B) < —0.25 mag, is 
P = 0.023 (2.0cr). 



7 A TOY MODEL FOR THE VIEWING 
ANGLE EFFECTS 

In the off-centre asymmetric SN la explosion scenario (§2), 
a correlation between Vnch and dm is expected. The viewing 

^ Here the dividing value, Ami5(B) = 1.13 is set merely by the 
requirement to split the sample into halves. 
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Table 4. Relations between Wneb ^^'^ d-m 
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Number of SNe. 

^ The best fit using the relation dm = a{v^^^/1000 km s~^) + /3. The errors for a and /3 are la. 

Probability that the distribution arises from a non-correlation. 

SNe with Ami5(B) ^ 1.12. 
e SNe with Ami5(S) 1.13. 
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Figure 10. The toy model used for the light curve calculations. 
The structure of the ECAP and HD zones is directly taken from 
the model of SN 2003hv (MlOa). 



angle effect on the peak brightness has been investigated 
by several authors. Sim et al. (2007) considered kinematic 
models for an off-centre distribution of ^^Ni, computed the 
radiation transfer in 2D, and showed that the luminosity 
can be different by ~ 0.5 — 1.5 mag for an observer at the 
offset direction as compared to the opposite direction. More 
recently, Kasen et al. (2009) investigated a series of off-centre 
hydrodynamic models in 2D, in which the deflagration was 
ignited in an off-centre manner (see also MlOc), obtaining 
peak magnitude differences at the level of ~ 0.5 mag for 
different viewing angles. There is also an expectation that 
the color is affected by the viewing angles, as investigated 
by Kasen et al. (2009) and Foley & Kasen (2010). However, 
there was no strong observational hint on the geometry of 
the inner ejecta at the time of their publications (but see 
Foley & Kasen 2010). In addition, they did not discuss the 
outcome of the off-centre ignition at the late phases as was 
investigated by MlOab. In the following, we investigate the 
effect of the viewing angle on the peak luminosity and color, 
adopting the simple explosion geometry proposed by MlOab. 
In the toy model shown in Fig. 10, the innermost regions 



are divided into two characteristic zones: a zone dominated 
by neutron-rich Fe-peak elements produced by electron cap- 
tures (ECAP; Electron-CAPture), the relatively high den- 
sity ^®Ni-rich region surrounding the ECAP region (HD; 
High-Density), and then the outer, relatively low-density, 
^^Ni-rich region on top of the HD region (LD; Low-Density). 
The ECAP and HD regions are suggested to be responsible 
for the formation of [Ni II] A7378 and [Fe II] A7155, respec- 
tively, from which we derived Uncb. Following MlOab, the 
ECAP and HD regions have offsets from the centre of 3, 500 
km s~^, whereas the LD region has only a slight offset in the 
opposite direction. The ECAP and HD zones are interpreted 
as products of the initial deflagration flame propagation, 
while the LD zone is a product of the subsequent detona- 
tion wave. There is possibly a variation of the offset velocity 
for different SNe la, although in the present analysis we use 
one specific value (3,500 km s^^) as our fiducial model|f| 
Note that this distribution was derived by modeling late- 
time spectra, but it is also qualitatively consistent with the 
outcome of hydrodynamic models if the defiagration sparks 
are ignited at an offset (MlOc). The hydrodynamic simula- 
tions have not yet explored all parameter space, nor been di- 
rectly compared to observations. We therefore take this 'ob- 
servationally constrained' toy model as our reference model. 

The masses and compositions of the ECAP and HD 
zones are taken from MlOa; the mass of is 0.01 Mq m 
the ECAP zone and 0.1 Mq in the HD zone. The mass of 
^^Ni in the LD zone is a parameter, corresponding to the di- 
versity of Mi^^m) and Ami5(B) for various SNe la. The ra- 
dial extent of the LD zone is also taken from MlOa. We calcu- 
late light curves for models with total M(^®Ni) from 0.21 to 
0.8 Mq, with 10% increases in M('^®Ni) between models. In 
the present analysis, we focus on the models with M (^^Ni) 
= 0.3 Mq (roughly consistent with the 0.3 — 0.4 Mq derived 
for SN 2003hv; Leloudas et al. 2009; MlOa) and 0.6 M© (for 
SNe la with a typical luminosity; see Stritzinger et al. 2006). 
We hereafter refer to the models as AO. 3, AO. 6, where the 
number denotes M(^''Ni). For all the models, the structure 
in the ECAP and HD zones (i.e., the masses of these zones 
and the offsets) was assumed to be exactly the same as the 
one derived for SN 2003hv. We fixed the masses in the ECAP 



* The presently available data analyzed by MlOab are consistent 
with no variation. To tackle the possible variation in the kinemat- 
ics and its relation to other explosion parameters would require 
a larger observational sample. 
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and HD zones since it has been suggested tiirougii late-time 
spectral modeling (Mazzali et al. 2007) that the masses of 
the innermost region are consistent with no variation for 
different SNe. 

The distribution of elements other than Fe-peak ele- 
ments was not discussed in MlOa, because a large fraction 
of them are located in the region above the LD-zone which is 
not accessible by late-time observations. These are assumed 
to be intermediate mass elements (IMEs) and/or a carbon- 
oxygen (CO) mixture. We assume an explosion energy of 
Ek = 1.4 X 10^^ erg and ejected mass of Mcj = 1.38 MqH 
These values determine the outermost velocity, as well as 
the density there. IMEs and CO are assumed to be dis- 
tributed uniformly in the LD and the outer zones. Accord- 
ingly, the mass in the LD zone is the sum of those of ^®Ni and 
IMEs/CO, whereas the outermost zone is purely composed 
of IMEs or CO. 

A question is whether or not the outermost region also 
shows an offset. To address this, we employ a simple model. 
It is assumed that the centre of the outer envelope is the 
same as for the LD-zone. This is motivated by the expec- 
tation that the outer envelope expansion is caused by the 
detonation wave. This configuration is roughly consistent 
with the finding by MlOb that SNe viewed from the offset 
direction (direction of the ECAP zone) are observed as LVG 
SNe while those from the opposite direction are HVG SNe 
because of the existence of a more extended envelope in the 
direction opposite to the offset. This characteristic is also 
seen in a hydrodynamic explosion model (see Ml Ob for a 
detailed discussion). 

For different masses of ^"Ni, we simulated the bolomet- 
ric light curve as a function of the viewing angle. The cal- 
culations were performed using the Monte-Carlo radiation 
transfer SAMURAI code. SAMURAI is a compilation of 3D 
codes that adopt Monte-Carlo methods to compute high- 
energy light curves and spectra (Maeda 2006a), optical bolo- 
metric light curves (Maeda et al. 2006b), and optical spectra 
from early (Tanaka et al. 2006, 2007) to late phases (Maeda 
et al. 2006c). For the light curve, the present version al- 
lows for a frequency-averaged gray opacity only. The optical 
transfer scheme follows prescriptions given by Lucy (2005) 
(see also Cappellaro et al. 1997). For the opacity to opti- 
cal photons, we adopted a phenomenological prescription of 
Mazzali et al. 

(20010 

Because of the simplified kinematic 
model and the limitations in the light curve calculations, the 
following results should be regarded as indicative. A more 



^ Strictly speaking, Ej^ is dependent on the final composition. If 
M(^^Ni) = 0.6 Mq and the remaining part is burned into inter- 
mediate mass elements, then Ek ~ 1-4 X 10^^ erg. By varying 
M(56Ni) between 0.2 and 0.8 Mq and assuming that the remain- 
ing fraction is burned into IMEs, the resulting Bk is in the range 
(1.3 — 1.5) X 10^^ erg. Thus, the variation in i?K is small, and we 
take -Ek = 1.4 x 10^^ erg as our fiducial value. 

We note that the rise time in the models is typically shorter 
than that estimated by observations. A part of the reason is likely 
the simplified opacity used in our calculations (see also Sim et al. 
2007). As wc are only interested in the 'relative' behavior, e.g., 
differences caused by different viewing angles and by different 
M(^^Ni), we did not try to obtain better agreement in the rise 
time. 
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Figure 11. The synthetic bolometric light curves for models AO. 3 
(blue-thick-solid lines) and AO. 6 (red-thin-solid lines). For each 
model, three curves are shown: Two angle-dependent light curves, 
with the angle measured from the direction of the offset of the 
ECAP/HD regions, and the mean angle-averaged light curve. For 
comparison, the mean light curves for models AO. 25 and AO. 4 arc 
shown (black-dashed lines). Time is defined as days past bolo- 
metric maximum. 



detailed and quantitative comparison between model and 
observations is beyond the scope of this paper. 

Figure 11 shows the bolometric light curves of mod- 
els AO. 3 and AO. 6, for different viewing angles. Model AO. 3 
shows a difference of ~ 30% in the peak luminosity depend- 
ing on the viewing angle. It is brightest if viewed from the 
offset direction, and faintest from the opposite direction. The 
reason for this is that photons from the HD zone have lower 
optical depth in the offset direction and preferentially escape 
into this direction. On the other hand, the contribution from 
the nearly spherical LD component is not very sensitive to 
the viewing angle. In model AO. 3, a large fraction of ^^Ni 
(~ 30%) is contained in the offset component. The bright- 
ness thus depends strongly on the viewing angle. On the 
other hand, in model AO. 6, the contribution of the offset 
component to M(^^Ni) is small {'^ 15%). Accordingly, the 
angle-dependence is small, at the ~ 10% level. 

Another important effect is that the light-curve shape 
can also be different for different viewing angles. Model AO. 3 
shows that the light curve is narrower and the rise time is 
shorter if viewed from the direction of the offset, while it is 
broader and the rise time is longer if viewed from the oppo- 
site direction (see also Sim et al. 2007). This stems from the 
shorter diffusion time scale in the direction of the offset. In- 
deed, this effect is part of the reason for the enhancement of 
the brightness in the offset direction, since faster rise to the 
peak results in a larger amount of radioactive 7-rays avail- 
able to power the peak brightness. The presently available 
data do not reveal this theoretically expected correlation 
between Ami5{B) and Vnch (MlOb), although the situation 
is still consistent with the idea presented here since a large 
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Figure 12. The behavior of the light curve shape and the peak 
luminosity of our models. The horizontal axis denotes Ami5 in 
the synthetic bolometric light curves, and the vertical axis the 
synthetic bolometric luminosity. Open circles denote angle aver- 
aged, mean light curves of different models (from AO. 21 to AO. 8). 
The angle-dependent behavior is shown for models AO. 3 and AO. 6 
by black filled symbols, where each point stands for the light curve 
from a certain direction. 



variation in Ami5(i3) is introduced by different amounts of 
^^Ni. With an increasing number of nebular spectra, this 
test will become possible. 

The trend in the peak brightness - light-curve width 
relation arising from different viewing angles works in a way 
opposite to the Phillips relation (see also Sim et al. 2007). 
A light curve that is brighter due to asymmetry is also nar- 
rower. For example. Fig. 11 shows that the peak brightness 
of the model AO. 3 viewed from the offset direction is com- 
parable to the mean (angle-averaged) brightness of model 
AO. 4, while the former results in a narrower light curve. 
Again, this effect is not as important for model AO. 6, for 
which the light curve shape does not strongly depend on 
the viewing direction. Figure 12 shows the peak luminosity 
in terms of Amis as derived from the synthetic bolometric 
light curves!"] 

Figure 13 shows the synthetic bolometric magnitude 
'difference' (~ dm) as a function of Wneb (i.e., the viewing 
angle) . It is compared to the V-band residuals as derived in 
the previous section. Given the uncertainty involved in the 
comparison (e.g., the simplified ejecta structure), it should 
be regarded as qualitative. 

The model behavior in Figure 13 was derived through 
the following procedures. First, we defined a bolometric 
Phillips relation for the models, based on the angle-averaged, 
mean light curves in Figure 12. Then, for a light curve of 
model AO. 3 and a given viewing direction, we searched for 



Note again that our simplified treatment in the opacity and 
in the ejecta geometry allows us to investigate only the relative 
behavior of different models (and viewing directions), but not to 
directly compare the models and the observational data. 



Figure 13. Model results compared to the observed sample in 
the J^nob vs. dm diagram. Model AO. 3 is shown by the shaded 
area, model AO. 6 by a solid line (see footnote 12). Points are for 
our sample of SNe, and a dashed line is the best linear fit to them 
(Fig. 8). 



a 'mean' light curve (spanning from model AO. 21 to model 
AO. 8) which shows a similar light curve shape. In general, 
the considered light curve lies between two mean curves of 
our model gridl^. Next, we define the model 'residual' as 
the magnitude difference between the angle- dependent lumi- 
nosity under consideration and the luminosity of the cor- 
responding mean light curve(s). This procedure is equiva- 
lent to deriving the observed residual (dm), but with the 
observed reference light curves replaced by the synthetic 
'mean' (angle-averaged) light curves. 

To compare the model residual (in bolometric magni- 
tude) and observed residual dm (in the V band), we have to 
make one assumption: The \^-band light curve is (at least 
qualitatively) traced by the bolometric light curve. Using 
a sample of SNe with the bolometric luminosity available 
(Contardo et al. 2000; Stritzinger & Leibundgut 2005), we 
confirmed that Knax is strongly correlated with the bolo- 
metric luminosity, and that the bolometric correction is not 
significantly different for different SNe. This justifies our as- 
sumption. 

In this comparison, not only can the dependence of the 
peak brightness on the viewing angle be important but also 
that of the light-curve shape. For instance, the residual for 
model AO. 3 viewed from the offset direction should be con- 
structed using the mean (angle-averaged) peak brightness of 
AO. 23 or AO. 25, which have similar light curve shapes, but 
smaller luminosity than, model AO. 3. A similar argument 

For model AO. 6, the light curve shape is not strongly depen- 
dent on the viewing direction (Fig. 11). Thus, for model AO. 6, 
the 'model residual' is computed as the magnitude difference be- 
tween the angle-dependent and mean light curves of the same 
model AO. 6. This is the reason why in Figure 13 the model resid- 
ual is expressed by a single line for model AO. 6, unlike model 
AO. 3. 
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applies to an observer in the direction opposite to the offset 
direction. Thus, the resulting model dm is generally larger 
than that estimated by simply taking the diflerence between 
the angle-dependent brightness and the mean brightness for 
the same model. For example, the residual for model AO. 3 
is ~ —0.4 mag if viewed from the offset direction, +0.3 
mag if viewed from the opposite direction. Thus the residual 
varies by up to 90 % (~ 0.7 mag) depending on the view- 
ing angle. This is much larger than simply taking the differ- 
ence between the angle-dependent brightness and the mean 
brightness for the same model, which is ~ 30% (Fig. 11). 

There is a qualitative agreement between the model 
and the data (Fig. 13). In particular, models with M(^®Ni) 
~ 0.3 Mq (model AO. 3) agree with the tentative correlation 
seen in the observed data. Also, if we calculate the expected 
difference in dm for SNe with blueshift (observed at < 90°) 
and redshift (> 90°), taking into account the solid angle, we 
obtain ~ 0.4 mag for this model. This is a bit larger than 
the observed difference in the two groups (Fig. 9 and re- 
lated discussion), but it explains the observed trend at least 
qualitatively. Note that the offset velocity has been directly 
constrained for model AO. 3 using SN 2003hv (MlOa). 

For brighter SNe la, model AO. 6 predicts a peak bright- 
ness almost independent of the viewing angle. We also find 
that the viewing angle effect on the light curve shape is not 
as important for this model. The reason for this is that most 
of the ^''Ni is in the 'spherical' LD zone in these models, since 
the amount of ^''Ni in the ECAP and HD zones (the "off- 
set" region) is fixed. Our analysis of the observational data 
suggests that the correlation is indeed weaker, as quantified 
by the correlation slope, for a subset of (brighter) SNe with 
smaller Ami5{B) than for those with larger Ami5(_B) (see 
Table 4 and related discussion). This is consistent with the 
model prediction, and may indicate that the contribution 
from the LD zone is more important and the viewing angle 
effect less important for brighter SNe. Although we fixed the 
mass of the offset component in our toy models, there is a 
theoretical expectation that the mass of the "offset" region 
is smaller for brighter SNe, since a larger asymmetry results 
in a weaker deflagration which is then followed by a stronger 
detonation in the off-centre deflagration-to-detonation tran- 
sition scenario (Kasen et al. 2009; MlOc). This will result in 
even smaller angle dependence for typical SNe la than model 
AO. 6 predicts. A larger sample will be required to clarify the 
relative contributions from the different zones, which in turn 
will provide hints on the geometry of SNe as a function of 
Ami5(B). 

Having computed angle-dependent luminosities, it is 
now possible to make a crude estimate on the relation be- 
tween the color and the viewing direction for our mod- 
els. The angle dependent effective temperature is approx- 
imately extracted from the model, by combining the lumi- 
nosity and photospheric velocity obtained by integrating the 
optical depth back inward along each radial direction (i.e., 
TeBF oc v'^^^L^^"^ where Wph and the luminosity L are an- 
gle dependent). Then we convert the effective temperature 
to the maximum color B — V assuming a relation between 
these two obtained for synthetic spectra by Nugent et al. 
(1995) (Their figure 3). Since the observational data (Fig. 
4) are already corrected for Amii,{B), we add a constant 
to the model B — V color (i.e., shifted vertically) so that it 
passes _B — 1/ ~ at I'ncb ~ (i.e., we focus on relative color 




Figure 14. Model results compared to the observed sample in 
the 

^neb VS. -Bmax — V^ax diagram (see the main text). Models 
AO. 3 and AO. 6 are shown by solid and dashed Unes, respectively. 
Points are for our low-extinction sample of SNe (Fig. 4a). 



change depending on the viewing direction). The result is 
shown in Figure 14. 

Since the model prediction will be sensitive to the as- 
sumed outer ejecta structure, this comparison should be re- 
garded to be qualitative. However, this does demonstrate 
that the basic features in the ejecta structure constrained 
by the observations (the inner structure from the late-time 
spectra and the outer structure from the relation between 
the early-time and late-time spectra; MlOab) predict the 
Vnch - color relation found in this paper as well, although 
the models have not been tuned at all to reproduce this re- 
lation. Indeed, Foley & Kasen (2010) have recently reported 
the same trend in the color and the viewing direction in their 
simulation based on a qualitatively similar ejecta structure, 
independently from this study. The qualitative behavior, i.e., 
bluer color for the direction closer to the offset, stems from 
the lower expansion velocity, thus higher temperature, in 
that direction. 

Unlike the bolometric luminosity, the relation is insensi- 
tive to the mass of ^^Ni, since the temperature is only weakly 
dependent on the luminosity. Under the assumptions within 
our toy model, the two models - AO. 3 and AO. 6 - do predict 
almost the same relation. This may indeed support why we 
see in the observations a strong relation in the color but not 
in the luminosity residual. 



8 DISCUSSION AND CONCLUSIONS 

In this paper, we have investigated how the explosion ge- 
ometry and viewing angle can influence the color and peak 
brightness of SNe la, and thereby lead to the residuals that 
remain in the peak magnitudes after application of the light- 
curve correction. We have used the wavelength shift of late- 
time emission lines to derive the viewing angle, and then 
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compared this quantity with the color and the luminosity 
residual for a sample of 20 SNe la. 



8.1 Intrinsic Color Variation 

We have found a correlation between the color at maximum 
and the wavelength shift seen in late-time emission lines. 
Using a sub-sample of 11 SNe which likely suffer from in- 
significant host reddening, selected based on the morpho- 
logical type of the host galaxy and the position of the 
SN within the host, we have found that SNe which show 
a blueshift/redshift in nebular emission lines (i.e., viewed 
from the offset / anti-offset direction in the off-centre ignition 
scenario) are bluer/redder than expected from the color- 
Amis (B) relation. This indicates that the previously sug- 
gested color difference between the 'high-velocity-gradient' 
(HVG) SNe and 'low-velocity-gradient' (LVG) SNe is con- 
nected to the viewing direction to the observer. 

In a next step, we expanded our sample, including ob- 
jects with potentially larger host reddening. Except for 3 
very red SNe [(Bmax- V;„ax) > 0.2 mag: SNe 1998bu, 2002bo 
and 2006X], all other objects agreed well with the color vs. 
Wneb relation derived from the "low-extinction" sample. This 
indicates that a significant part of the color excess previously 
attributed to host extinction is actually due to intrinsic color 
variations. 

This raises the question whether the often preferred 
value of Rv < 2 (much smaller than the typical Galactic 
value of ~ 3.1) really refiects different properties of inter- 
stellar/circumstellar dust around SNe la. FIO argued that 
the optical color minus NIR color relation can be reproduced 
even with Rv ~ 3.2 once heavily reddened SNe are omitted. 
On the other hand, they required Rv ~ 1 — 2 to minimize the 
dispersion in the standardized luminosity calibration. They 
pointed out that an intrinsic color variation which does not 
correlate with the decline rate, but does correlate with the 
luminosity, may solve this apparent discrepancy. Indeed, the 
intrinsic variation related to the viewing angle that we have 
found in this paper does not correlate with the decline rate. 
If we exclude SNe which are clearly heavily reddened, we 
see that Rv as large as the Galactic value could be accept- 
able. This issue requires more careful and systematic analy- 
sis based on a larger sample. Also, we note that even if this 
works for the majority of SNe la, some heavily reddened SNe 
appear to require Rv ~ 1 — 2, as highlighted by SN 2006X 
(FIO). 

A toy model constructed with the constraints so that it 
explains the late-time spectral variation (MlQa) and its re- 
lation to the early-phase spectral diversity (MlOb) predicts 
the bluer/redder color for smaller/larger (blueshift/redshift) 
Wneb, as is consisteut with the observational data. For our 
model sequence, the predicted trend is not sensitive to 
M{^^m) (or Ami5(B)). 

8.2 Second Parameter in the Luminosity 
Calibration? 

We have also investigated a relation between the SN la lu- 
minosity residual after calibration with the Phillips rela- 
tion, and the wavelength shift seen in the late-time nebu- 
lar emission lines. For our small sample, the correlation is 



not strong, and we regard this result as tentative. Keeping 
this caveat in mind, we see a tendency that SNe la with a 
blueshift/redshift in the nebular lines show a higher/lower 
peak luminosity than expected by the Phillips relation. 
There is an average difference of ~ 0.25 mag in the peak 
magnitudes between SNe la showing a blueshift and those 
showing a redshift. 

Calculating the light curves based on the geometry de- 
rived by MlOa for the relatively faint SN la 2003hv, we have 
found that SNe la viewed from the offset direction should 
have a peak brightness larger than the mean, and those 
viewed from the opposite direction should be fainter. The 
difference between these two extreme observer directions is 
~ 0.7 mag in our fiducial model AO. 3B This behavior is 
consistent with the observational data. Also, the averaged 
difference between SNe showing blueshifts and redshifts is 
~ 0.4 mag in this model, enough to explain the observed 
value (~ 0.25 mag). 

The comparison between SNe la with normal/large 
peak luminosity and the models is less straightforward, since 
the geometry of such explosions has not been directly con- 
structed from observations. Assuming that the degree of 
the off'set is similar to that of SN 2003hv, we expect that 
the viewing angle effect is less pronounced for brighter SNe 
la. This is consistent with the behavior we have found in 
the data. Further variation is expected since the degree of 
the ofi^set could be difi^erent for SNe with different /S.miz,{B) 
(e.g., Kasen et al. 2009; MlOc). We therefore suggest that 
the residual could be explained by a combination of the 
configuration [e.g., the relative contribution between the 
ECAP/HD and LD zones, which may well be expressed by 
one parameter, i.e., Ami5(i3)] and the viewing angle. Thus, 
we do not expect a single straight line to give a perfect fit 
to the residual vs. late-time velocity plot. This may be one 
reason why the correlation in Fig. 8 is not too strong. 

8.3 Future Perspectives and Cosmological 
Applications 

Our results could be further tested by polarization measure- 
ments. SNe la generally show small continuum polarization 
around maximum brightness, indicating that they are more 
or less spherical (e.g., Wang et al. 1996). However, early- 
phase polarization measurements mainly probe a region near 
the surface of the SN la ejecta, which we also suggest to be 
nearly spherical. Therefore, the small polarization is likely 
not a strong argument against our present interpretation. 
The low continuum polarization is sometimes translated into 
a small deviation of the photosphere from spherical symme- 
try, and hence a small dependence of the brightness on the 
viewing angle (e.g., Hoflich et al. 2010). However, this state- 
ment depends strongly on the assumed geometry. For exam- 
ple, a continuum polarization of < 0.3%, as is usually found 
in SNe la, implies an axis ratio of less than 10% for an ellip- 
soidal photosphere (Hofiich 1991). However, for a one-sided 
distribution of ^®Ni, as suggested in the present work, the 
expected continuum polarization is generally much smaller 
than for an ellipsoid (Kasen & Plewa 2007) despite a larger 

^■^ Including the effect that the light curve shape is different for 
different viewing directions. 
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expected variation in the angle-dependent brightness than 
for the ellipsoidal case (e.g., Sim et al. 2007; Kasen et al. 
2009; this work). 

There is a correlation between the velocity gradient and 
the Si II line polarization (Leonard et al. 2005; Chornock & 
Filippenko 2008; Maund et al. 2010). This may indicate that 
the viewing direction is indeed controlling the line polariza- 
tion level, if the velocity gradient is determined by the view- 
ing angle effect (MlOb). Further study on the polarization 
of SNe la both in observations (e.g., Wang et al. 2007) and 
in theory (e.g., Hoflich 1991; Kasen & Plewa 2007) should 
provide a good test for the geometry of SNe la, and thus for 
the results in the present work. 

A problem in our analysis, especially for the residual 
issue, is the small sample size. The uncertainty in the dis- 
tance measurement is also critical in our investigation of the 
viewing-angle effect on the luminosity residual. We suggest 
the following strategies to decrease the uncertainty in the 
distance estimate; (1) obtain late-time spectra for SNe la 
at redshift z > 0.02. At z 0.02, the V-band magnitude 
of typical SNe la at ~ 200 days after maximum bright- 
ness is ~ 21 — 22 mag (depending on the extinction). Spec- 
troscopy is thus possible with 6 — 8m class telescopes, and 
(2) obtain comprehensive photometry during both the peak 
and tail phases. The residual arising from the viewing angle 
can in principle be seen in the peak-to-tail luminosity ratio, 
since the effect of the viewing angle vanishes at late-phases 
(e.g., Maeda et al. 2006b). This may provide a distance- 
independent measurement of the residuals caused by the 
viewing angle effect. Such observations are less demanding 
than spectroscopy, and thus can reach to larger redshift. 

One of the current limitations in SN la cosmology is 
the fact that a dispersion at the level of ~ 0.15 mag remains 
after the standardization of the peak luminosity with exist- 
ing relations, and that the physical origin of the residual has 
not been identified. The viewing angle effect may account for 
part of the dispersion of the SN la luminosity calibration, 
although the presently available sample does not allow us to 
quantify how much improvement can be achieved by taking 
this effect into account. However, the effect of the random 
viewing angle enters in the SN la luminosity calibration as 
a source of a statistic error. Thus, increasing the number of 
SNe la for cosmology should effectively reduce this error in 
estimating cosmological parameters. 

In addition, it may be worthwhile looking into the fre- 
quency distribution of the residuals. Although we expect 
that the statistical error is introduced by random viewing 
angles, there is no reason to expect it to obey a Gaussian 
distribution. Investigating a non-Gaussian component in the 
scatter of the Hubble diagram may provide additional in- 
sight. For example, a larger SN la sample may allow us to 
estimate the non-Gaussian contribution to the statistical er- 
ror, which will then provide a quantitative estimate of the 
viewing angle effect independent of uncertainties in the dis- 
tance measurements for individual SNe. 

The result of the present work may shed light on how 
to develop a more accurate SN la cosmology than currently 
employed. The correlation between the color near maximum 
and Vnch enables us to discriminate the intrinsic color and 
the host extinction when late-time spectra are available. 
This opens up a possibility of studying the intrinsic color 
and the host extinction, including Rv, in detail for nearby 



SNe up to a redshift of ~ 0.02 for which late-time spec- 
troscopy is possible. This will hopefully provide information 
about how to distinguish the intrinsic color and host extinc- 
tion, the information applicable to high-redshift SNe. The 
relation between the velocity gradient and the nebular emis- 
sion line shift (MlOb) suggests that one could use the veloc- 
ity gradient (accessible to higher redshift SNe) instead of 
the latter, in the color and luminosity calibrations (see also 
Foley & Kasen 2010). Indeed, it has been argued that LVG 
and HVG SNe should be treated differently in the luminosity 
calibration (e.g., Wang et al. 2009b). Further investigating 
relations among these observable (e.g., Fig. 5), as well as 
finding other observables which correlate with i^ncb, could 
be useful in improving the color and luminosity estimates 
also for high- 2: SNe la. 
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data provided by Rieke & Lebofsky (1985)]. For SN 2001el, 
the value of Av was taken from Krisciunas et al. (2007) who 
used SN 2004S as a template because of the similarity be- 
tween these two events. When these near-IR measurements 
were available, we adopted Ay based on these estimates. 
When near-IR measurements were not available, we con- 
verted E{B — V) color excess estimates from the literature to 
Av using Rv ~ 1.72. In the cases of SNe 2007sr (Schweizer 
et al. 2008), 2007on and 2009ab (Stritzinger et al. in prep.), 
E{B — V) was obtained using the template light curve fitter 
SNooPy (Burns et al. 2011). The values thus compiled are 
listed in Table Al. 

From this comparison, we find the similar results to 
those obtained by the comparison to E{B — V) derived by 
Wang et al. (2009b). Our estimates tend to be smaller, and 
even if we assume Rv ~ 3 to convert our E{B — V) estimate 
to Av, the extinctions we derive are mostly consistent with 
those derived by the other method (except for the heavily 
reddened SN 2006X). 



Figure Al. Comparison of the extinction. The horizontal axis 
denotes the values of Av derived in this study (Tab. 2), while 
the vertical axis provides values from literatures (Tab. Al). The 
symbols and color coding indicate the velocity gradient near max- 
imum brightness (see Fig. 4 caption). 
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APPENDIX A: EXTINCTIONS AS COMPARED 
TO OTHER ESTIMATES 

In §4, we compared the color excesses we derived (with 
Vncb) to those estimated by Wang et al. (2009b). As a fur- 
ther, additional test, we also check how our extinction val- 
ues compare to those previously estimated and presented in 
the literature by different methods. Assuming Rv = 1.72 
(§1, §5), we convert E{B — V) to Av- Figure Al shows a 
comparison between Ay as derived in this study and the 
values from the literature. The literature values were se- 
lected as follows: if optical minus near-IR (i.e. V — K, V — J, 
V — H) color curves are available, the template color curves 
of Krisciunas et al. (2000, 2001, 2009) were used to de- 
rive Av with the infrared extinction law given by Rieke & 
Lebofsky (1985). A similar estimate is also possible using 
the near-IR spectral energy distribution (SED) at a single 
epoch [e.g., Wang et al. (2008), who used the extinction law 
of Cardelli, Clayton, & Mathis (1989) constructed from the 



APPENDIX B: UNCERTAINTIES IN THE 
LUMINOSITY CALIBRATION 

Compared to the color calibration (§4), our analysis on the 
luminosity residual (§6) suffers from various sources of un- 
certainties. These are independent from i^ncb, and so it is 
unlikely that any correlation we investigate in the present 
study is artificially introduced by these uncertainties in a 
systematic way. However, our present sample is still small 
and thus statistical errors due to these uncertainties can still 
be important. In this section, we investigate how our results 
in §6 are affected by changing the procedures to estimate 
the distance, reddening, and standardized luminosity (§5). 
Results are summarized in Table Bl. Each item in the table 
is explained in the following. 

Bl Distances 

To examine the uncertainty in the distances, we first com- 
pare the distances obtained by different methods, when 
available, in Fig. Bl. Hubble fiow (HF) distances with a re- 
cession velocity below 1, 000 km s^^ are not considered as al- 
ternatives here. Different measurements typically agree with 
each other to within the errors (albeit these errors are some- 
times quite large, when only the HF distances are available 
as the 'other' measurement). However, there are several ex- 
ceptions. NGC 4639 (SN 1990N) and NGC 4321 (SN 2006X) 
are (possible) members of the Virgo cluster. If we adopt the 
distance modulus ^ = 31.17±0.04 mag for the Virgo cluster 
(Kelson et al. 2000), it differs significantly from the Cepheid 
distance to NGC 4639 we adopted {^i (KP) = 31.71 ± 0.15 
mag). We note, however, that Riess et al. (2009) revised the 
Cepheid distance to NGC 4639 to be ~ 31.48 ± 0.13 mag. 
Although there is still a discrepancy between the distances 
to NGC 4639 and the Virgo cluster, this new measurement 
makes the agreement better. If we adopt the value of Riess 
et al. (2009), SN 1990N should be fainter by 0.23 mag than 
our fiducial estimate. We did not adopt this value in the 
main text, to avoid possible systematic errors in different 
Cepheid measurements. The distance to NGC 524, the host 
of SN 2000CX is 31.74 ±0.20 mag (SBF) or 32.60 ±0.27 mag 
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Table Al. Extinction from the literature 



SN 


Ay 

(mag) 


Method 


Reierences 


1990N 


0.16 




Ph99 


1994D 


0.0 




Ph99 


1997bp 


0.46 




A04 


1998aq 


0.0 




W09b 


1998bu 


1.06 


NIR 


KOO 


2000cx 


0.0 




LOl 


2001cl 


0.47 


NIR 


K07 


2002bo 


0.51 


NIR 


K04 


2002dj 


0.0 


NIR 


P08 


2002or 


0.40 




W09b 


2003du 


0.0 


NIR 


S07 


2003hv 


0.0 




L09 


2004dt 


0.19 




W09b 


2004eo 


0.0 


NIR 


P07 


2005cf 


0.18 


NIR 


W09a 


2006X 


2.1 


NIR 


W08 


2006dd 


0.12 




SlOa 


2007on 


0.25 




SlOb 


2007sr 


0.13 




SlOb 


2009ab 


0.076 




SlOb 



" Refs. A04. Altavilla et al. (2004); KOO. Krisciunas et al. (2000); K04. Krisciunas et al. (2004); K07. 
Krisciunas et al. (2007); LOl. Li ot al. (2001); L09. Leioudas et al. (2009); P07. Pastorello et al. (2007a); 
P08. Pignata et al. (2008); Ph99. Phillips et al. (1999); S07. Stanishev et al. (2007); SlOa. Stritzinger et 
al. (2010); SlOb. Stritzinger et al. (in prep.); W08. Wang et al. (2008); W09a. Wang et al. (2009a); W09b. 
Wang ct al. (2009b) 



Table Bl. Relations between t^ncb and dm. 



Description 






oL 






/3 


pc 


significance 


Distance'^ 


20 


0.092 : 


± 0..052 


0.070 


± 


0.077 


0.038 


1.80- 


STS^ ~ 
RvL 


20 


0.11 


± 0.050 


-0.013 


± 


0.076 


0.016 


2.10- 


20 


0.072 


±0.055 


0.040 


± 


0.082 


0.099 


1.3(T 


Ay (literature^ 


20 


0.050 


± 0.049 


-0.064 


± 


0.074 


0.16 


l.Oo- 


Mv(Ami5(S):ii 


20 


0.081 


±0.051 


0.071 


± 


0.077 


0.059 


1.6ct 



" The number of SNe. 

^ The best fit using the relation dm = a (Dn(,b)/1000 km s"-*^) + /3 mag. The errors for a and /3 are Icr 
uncertainties. 

Probability that the distribution arises from a non-correlation. 
<^ Obtained by changing the distances to SNe 1990N and 2006dd (§B1). 
^ Obtained by changing the distances to the STS scale (§B1; Table B2). 

/ With Rv = 3.1 (§B2). 

S' Ay adopted from the literature (§B2). 

^ The reference magnitude changed to that of Phillips et al. (1999) (§B3). 



(HF), which clearly do not agree with each other. Although 
we usually adopt the SBF distance as the better estimate, 
we note that another supernova, SN 2008Q, appeared in the 
same galaxy and both of these SNe would be peculiar out- 
liers in this case. Thus, we adopted the HF distance to SN 
2000CX (§3). Finally, the PNLF distance (/x = 31.26 ± 0.1 
mag) to NGC 1316 (SN 2006dd) is smaller than the SBF 
distance (p = 31.50 ± 0.17), and Stritzinger et al. (2010) 
argued that the former is more accurate (thus adopted in 
the main text). 



Following the above inspection, we change the distance 
to SN 1990N (from the original KP value to that derived by 
Riess et al. 2009), and that to SN 2006dd (from the PNLF to 
the SBF), and repeat the same analysis as we did in Fig. 8. 
We thereby obtain a chance probability of P = 0.038 (1.8(t) 
(Table Bl). 

As another test, we explore how our results are affected 
if we adopt the STS Cepheid distances (§5). For 4 SNe la 
with a Cepheid distance available, there are also measure- 
ments by the STS group (Saha et al. 2006). We also change 
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Figure Bl. Comparison of the distance we adopted (fio) and 
that from a different method (/^other) for SNe with several inde- 
pendent distance measurements available. SNe discussed in §B1 
are marked by red-filled circles and indicated by SN (/iq — /Mother)- 



Figure B2. The same as Fig. 8, but using the STS measurements 
for the Cepheid distances and a revised zero point for the SBF 
distances. See the caption of Fig. 4 for the meaning of the different 
symbols. The dotted line is the best-fit line for the STS distances, 
while the dashed line is that for the KP distances. 



Table B2. SNe la sample on the STS distance scale. 



SN 



Ifost 



(mag) 



1990N NGC 4639 32.20 ± 0.09 (STS) 

1994D NGC 4526 31.18 ± 0.20 (SBF) 

1998aq NGC 3982 31.87 ± 0.15 (STS) 

1998bu NGC 3368 30.34 ± 0.11 (STS) 

2002bo NGC 3190 31.90 ± 0.24 (SBF) 

2003hv NGC 1201 31.57 ± 0.3 (SBF) 

2006X NGC 4321 31.18 ± 0.05 (STS) 

2007on NGC 1404 31.65 ± 0.19 (SBF) 



the SBF measurements to be consistent with the STS mea- 
surements. The SBF values used in the main text are cah- 
brated with the KP Cepheid distance zero-point. Since the 
STS values are on average larger than the KP values by 
~ 0.2 mag, we add 0.2 mag for the SBF distances. Table B2 
lists the STS distances thus derived for SNe la. For SNe la 
missing in Table B2, the same distances as in Table 1 are 
used. Figure B2 shows the ^ncb vs. dm diagram with the STS 
distance scale. For 20 SNe la, we obtained a chance probabil- 
ity of P = 0.016 (2.1(j). The slope of the linear fit is steeper 
than that derived with the KP distance scale, although they 
are consistent within the errors of the fits. The zero-point 
is displaced by ~ 0.1 mag, reflecting the difference between 
the STS and KP distance scales. 



B2 Extinction 

Another large uncertainty comes from the estimate of the ex- 
tinction within the host galaxies. To check the uncertainty, 
we change the value of Ry to 3.1, the standard value for 
Galactic extinction. Following §4 and FIO, Rv for SN 2006X 



is left unchanged. While it has been argued that the typi- 
cal Galactic value is not necessarily applicable for the host 
galaxies of SNe la (e.g., Folatelli et al. 2010; Hicken et al. 
2009ab; Wang et al. 2009b; Yasuda & Fukugita 2010), our 
analysis in §4 suggests that an Ry close to the typical Galac- 
tic value can be acceptable for mildly extinguished SNe, once 
the intrinsic color variation due to the viewing direction is 
taken into account. Figure B3 shows the v^eb — dm diagram 
with Rv = 3.1. We obtain linear fit parameters as follows; 
P = 0.099 (1.3cr). The best fit lines for different Rv (1.72 
and 3.1) are consistent to each other within the fitting errors. 
The significance of the fit is slightly weaker for Rv = 3.1 
than in our fiducial case. It is, however, a reasonably large 
range of Rv that we investigate. 

Next, we explore whether our use of the relation be- 
tween Vnch and the color to derive the extinction might affect 
the results. For this purpose, we replace the host galaxy ex- 
tinction by the values given in the literature, and repeat the 
same analysis. For the literature values, see Table Al and 
the related discussion in the main text. The result is shown 
in Figure B4. The fitting resuff is P = 0.16 (l.Oa) (Table 
Bl). Although there is still a correlation, it is weaker than 
if the Vncb " color relation derived in this paper is adopted 
(1.6a). 



B3 Standardized Luminosities 

By definition, dm depends on the standardized luminos- 
ity, i.e., Mv(Ami5(B)). We have used an updated relation 
given by FIO. To check the related uncertainty, we replace 
My (Amis (P)) by the original version of the Phillips rela- 
tion including the second order term in Ami5(P) (Phillips et 
al. 1999) and repeat the same analysis. My (Ami5(P) ~ 1-1) 
is set to be —19.12 mag, to be consistent with Folatelli et al. 
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Figure B3. The same as Fig. 8, but using Ry = 3.1 for the ex- 
tinction. See the caption of Fig. 4 for the meaning of the different 
symbols. The dotted line is the best-fit line for the data points 
with Rv =3.1, while the dashed line is for the fiducial value of 
Rv = 1.72. 



Figure B5. The same as Fig. 8, but using the original Phillips 
relation (Phillips et al. 1999) for the reference magnitude 
(A?v(Ami5(_B))). See the caption of Fig. 4 for the meaning of 
the different symbols. The dotted line is the best-fit line for this 
data set. 
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Figure B4. The same as Fig. 8, but Ay replaced by the values in 
the literature (Table Al). The best-fit line is shown by the dotted 
line. 

(2010). Figure B5 shows the resulting v-nch — dm diagram. 
We obtain P = 0.059 (1.6cr). 



formed a linear regression to a set of variables (Xiiaxi , Yi ± 
(Ty. ) (where i spans the numbers in the SN sample). First, we 
produced 10^ test distributions obtaining (X,' j.,y/j.) (with 
k spanning from 1 to 10^), where a Gaussian distribution is 
assumed for the variation in Xi and Yi with the associated 
errors ax^ and ay. . Here, cry. includes both the extinction 
and distance uncertainties. For each fc-th test distribution, 
we performed a linear regression fitting assuming the func- 
tional form Y' — UkX' + Pk- Here we used only the error 
associated with the extinction as a weight in the fitting 
for each distribution, because the 'relative' errors between 
different distance measurements are difficult to quantify. We 
thereby obtained at = Ofe.o ± (^a^ and fit = /?fe,o ± o-0^ by 
the linear regression for each fc-th test distribution. This re- 
sults in a probability distribution of the fitting coefficients, 
a and j3, by convolving 10^ Gaussian distributions for Ok 
and Pk- The final fitting result is then obtained by fitting 
Gaussian profiles to the probability distribution of a and /3. 
To estimate the chance probability P that the distribution 
{X ± ax, Y ± cry) would arise from a non-correlation, we 
counted the probability for a ^5 (if the mean value of a is 
positive). 

This paper has been typeset from a TJ^X/ M^jX file prepared 
by the author. 



APPENDIX C: MONTE-CARLO SIMULATION 
FOR ESTIMATING THE CHANCE 
PROBABILITY 

To estimate the chance probability P that a distribution 
arises from a non-correlation, we have used the following 
method based on Monte-Carlo simulations in §6.: we per- 



